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Bound states of heavy qq quarks are reviewed within the context of QCD, paying attention to 
what can be derived from the theory with a reasonable degree of rigour. This is compared with the results 
of semiclassical arguments. Among new results, we report a very precise 0(oQ) evaluation of b, c quark 
masses from quarkonium spectrum with a potential to two loops. 



1. INTRODUCTION 

In the present note we are going to review some as- 
pects of the QCD analysis of heavy quarkonia, cc 
and especially bb states. This is fitting for a con- 
ference which (slightly ahead of time) celebrates 
the 25th anniversary of QCD. Indeed, the theory 
of quark interactions became a respectable theory, 
QCD, only with the advent of asymptotic free- 
dom in 1973. Before that date we had the quark 
model, a somewhat inconsistent set of semiphe- 
nomenological calculations. Among these an im- 
portant role was played by bound state calcula- 
tions in the so-called constituent quark model, 
developed in the early sixties by, among others, 
Morpurgo, Dalitz and collaborators, and Oliver, 
Pene, Reynal and Le Yaouanc. In this model 
u,d, s quarks were given phenomenological masses 

of 300 500 MeV, and were bound by potentials: 

the harmonic oscillator potential being a popular 
choice because of its simplicity. Quite surprisingly, 
a large number of properties of hadrons could be 
reproduced in this way. 

After the advent of asymptotic freedom, 
and with it a consistent field theory of strong in- 
teractions, it was possible to show that, at least 
for heavy quarks and at short distances, the inter- 
action is of Coulombic type. For colour singlet qq 
states of the form 

(1.1) 

r 

In one of the first applications of QCD, De Rujula, 
Georgi and GlashowM showed that taking into 
account relativistic corrections and colour alge- 
bra one could calculate the spectrum of the then 



known hadrons, including in particular such fea- 
tures as the N — A splitting, and even the E° — A 
splitting, something that had defied previous, non- 
QCD analyses. They were also able to predict 
some qualitative features of the charmonium spec- 
trum. 

Nowadays we expect more from QCD, at 
least for heavy quarks. The reason is that there it 
can be easily proved that, to leading order in (v 2 ) 
(with v the velocity of the quarks) the interaction 
can be described by a potential. At very short 
distances this potential has to be of the Coulombic 
type, Eq. (1.1); but even at long distances the 
corrections to this are expected to be of the form of 
a function U(r). At short distances (1.1) should be 
modified by radiative corrections, but these should 
be of the form of a function of r. 

Needless to say, relativistic corrections will 
in general not be representable by potentials, as 
is the case even in QED. In QCD one encounters 
QED-like corrections and idiosincratic QCD ones 
associated with the complicated structure of the 
vacuum. In particular we have those involving the 
gluon condensate (a s : G 2 :), first studied in this 
context by Leutwyler and VoloshW 2 ! (the quark 
condensate also gives contributions, but, for heavy 
quarkonium, subleading ones). 

2. SHORT DISTANCE QUARKO- 
NIUM: PURE QCD ANALYSIS, b AND 
c QUARK MASSES 

For very heavy qq bound states the equivalent of 
the Bohr radius, a = 2/(mCi7Q s ), is much smaller 
than the confinement radius, R ~ A" 1 . So we ex- 
pect that, for lowest n states, with n the principal 
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quantum number, confinement may be neglected, 
or at least treated as a first order perturbation. 
In this case the potential may be obtained from 
perturbative QCD. At tree level (Fig. 1) we get 
the Coulombic potential, Eq. (1.1). Including ra- 
diative corrections will yield improved approxima- 
tions, in particular giving a meaning to the quan- 
tity a s in (1.1). 




FIGURE 1. One- gluon ex- 
change. 

These radiative corrections have been eval- 
uated by a number of people. Those to the spin- 
independent part of the potential, in the strict 
static approximation, were first calculated by Bi- 
lloire' 3 '. Relativistic corrections were evaluated 
in refs. 4, 5 and they were completed in ref. 6 
where also some pieces of the two-loop correc- 
tions were given. A partial evaluation of the 
static two-loop interaction has been published by 
Peter while the completed calculation has been 
performed very recently by the author ^ for the 
n = 1, I = state, and will be given here for the 
first time. 




FIGURE 2. Some radiative correc- 
tions. 

To take into account all terms giving cor- 
rections of 0{a 2 s ) to the energy spectrum one 
writes the Hamiltonian as 

H = H {0) +H 1 (2.1a) 

where H^ may be solved exactly and contains 



the Coulomb-like part of the interaction: 

g(°>=2m+^A- C ^ 2) 



7eA)\ a s (M 2 ) 



+ 



7e I ai(3o 



Pi 



l2 +7E ~ + bl 



7T 

a s 



(2.1b) 



H x is 

Hi = \4rcc + vl L) + V} L) + V (LL) + V S . IC 1 + V hi , 

(2.1c) 



-1 2 C F a s 

nroo — -. jfl H o — 

4r7j J m z r 



A, 



-C F o?. 



-C F a s (^, 2 ) 2 0o logr^ 



V, 



(L) 



aiPo + -T- H 7T 



y{LL) 

V h{ = 
In above equations, 



Cf(3o^ s l°g r M 



47T 2 
Ks.rel 

AixC F a, 



r 

C F a 2 a 2 s 



3m 2 



2mr 2 ' 
«(s + l)<f(r). 



31C A -20T F n f 
a i = — ~ — 1-47; 



36 

C F — 2Ca 



~ -2.33; 



bi 



16 { [' 16:: 1 " : " :i 

_ [igt + f C(3) ] C A T F n f - 
[f _16C(3)]C7 F T F n / + flr 2 n 2 }, 



with a\ calculated in ref. 3, b\ in ref. 7 and ai and 
many of the rest of the terms in ref. 6. All terms in 
Hi are to be treated as first order perturbations 
of except for the term , which has to 

be treated to second order. Thus it produces, in 
addition to the first order contribution, 



m 4tt V ° S 2 



7e 



(2.2) 
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the second-order energy shift. For the ground 
state it is very small, of about 4 MeV. 

The first order contributions of the other 
V's are easily evaluated using the formulas of 
ref. 6. 



R 




FIGURE 3. The region where the 
quark pair move inside the con- 
finement region. 

The leading nonperturbative (NP) correc- 
tions can be shown to be those associated with 
the contribution of the gluon condensate. They 
may be understood as follows. We consider that 
the quarks move in a medium, the QCD vacuum, 
which is full of soft gluons (Fig. 3) that we rep- 
resent by their field strength operators, G° (x). 
When a <C R, we may consider that the confine- 
ment size is infinite and, moreover, one can ne- 
glect the fluctuations of the G° (x) in the region 
of size a in which the quarks move. So we may 
approximate the effect by introducing an interac- 
tion, which in the static limit will be of dipole 
type, of the quarks with a constant gluonic field, 
H NP = t c r t G c oi (0)- We consider that (G^) = 0, 
but (a s : G 2 :) ^ 0. For dimensional reasons, 
this will give the leading NP contribution to the 
spin-independent energy shifts, which are of the 
formal , 

7ro 6 e n t(q a : G 2 :) 
6spE,a - m {mCFas)i , (2-3) 



where the numbers e n i are of order unity, eio — 
1.5. The evaluations for the spin-dependent shifts 
may be found in the second paper of ref 8 (with a 
correction in ref. 9) and the contributions of higher 
order operators has been considered in ref. 10. 
Note that, as already remarked by Lcutwylerl 2 !, 



one cannot derive (2.1) from a local potential; but 
the effect may be approximated by a cubic one, 

^Gluoncond.W ~ A 3 . (2.4) 

Let us ummarize the results! 6 ' 8 ' 9 ' 10 !. The 
calculation is fully justified, in the sense that 
higher order corrections (both perturbative and 
NP) are smaller than lower order ones for bb with 
n = 1. The same is partially true for the energy 
levels of the same states with n = 2 and, for cc, 
for n = 1. For the wave functions of bb, n > 2 
and all cc states, and for the energy levels with 
higher values of n than the ones reported above, 
the calculation is meaningless as nominally sub- 
leading corrections overwhelm nominally leading 
ones. 

For bb one gets a precise determination 
of nib and rh b (m 2 ) (pole and MS masses), a re- 
liable prediction for the hyperfine splitting, and 
reasonable agreement with the experimental value 
of Y — > e + e~: 

m b = 4906+12 MeV 

T 1S (2.5a) 
fh b (fhl) = 4397t^ MeV, 

r{T -> e+e - ) = 1.12 keV (cxp : 1.32 ± 0.05). 

(2.5b) 

For cc a reasonably accurate value for is also ob- 
tained for m c : not including the estimated sys- 
tematic error, 



m c = 1570 ± 20 MeV 
m c {fh 2 c ) = 1306^ MeV. 



(2.6) 



These results are obtained with the one-loop po- 
tential with relativistic corrections' 6 '. We may ex- 
tend the calculation to two loops, using the Hamil- 
tonian of Eqs. (2.1,2) above, plus leading NP cor- 
rections, Eq. (2.3). Taking A = 200 MeV, the 
renormalization point /x = 2/a ~ 2.5 GeV, and 
varying /i 2 by a factor two to get the systematic 
errors of the calculation one finds from the T and 
J /if) masses the (pole) quark masses ^ correct up 
to, and including, 0{a A s ) terms: 

m b = 4984 ± 62 MeV, m c = 1797 ± 70 MeV . 

— ^ 
The corresponding MS bar masses are m,b(m 2 ) = 

4.446 GeV and m c (m 2 ) = 1.501 GeV. The val- 
ues of the masses are slightly larger than those 
one finds with the sum rule method (see for ex- 
ample, refs. 12). This may be easily understood 
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if one realizes that the last are obtained in calcu- 
lations accurate to 0(a 2 ) while the ones reported 
here include terms in o? s (Eq. (2.5)) and a*, for 
Eq. (2.6). If we had only included the terms in a 2 
in a potential calculation we would have obtained 
m& = 4746 MeV, for example. This is comparable 
to the sum rule value, so the discrepancy is seen to 
lie in the contribution of terms of order a 3 , not 
taken into account in the sum rule evaluations. 

3. QUARKONIA AT LONG DIS- 
TANCES. CONNECTION BETWEEN 
THE LONG AND SHORT DISTANCE 
REGIMES 

Here we consider bound states of heavy quarks 
at long distances. This certainly includes cc with 
n > 1 and bb with n > 2; n = 1 for the first 
and n = 2 (and, a fortiori, n = 1) for the sec- 
ond are somewhat marginal. As stated in the 
previous section, perturbative QCD supplemented 
with leading NP effects fails now; but, fortunately, 
and since the average velocity of bound states de- 
creases with increasing n, we expect the dynamics 
to be governed by a potential: our task is to de- 
termine it. This has been considered by a number 
of people! 13-17 ! Here we will follow the derivation 
of ref. 16 in the version of rcf. 18, wich will allow 
us to establish connection with the short distance 
analysis of the previous section. 

The potential, that we denote by V(r), is 
expected to exhibit a number of features. First 
of all, it should behave as or at long distances. 
Secondly, it should contain a Coulombic piece, so 
we write 

V{r) = -Z + U(r), (3.1) 

and, at short distances, one should be able to iden- 
tify k = Cfol s + radiative corrections. 

To find this potential consider the Green's 
function in terms of the Wilson loop, working di- 
rectly in the nonrelativistic approximation, and 
for large time T: for a qq pair: 

G(x,x;y,y) = J VzVze-^ Ko+R °\W(C)), 

(3.2) 

with Ko,Ko the kinetic energies, 

K = ^ / dz(t) 2 , etc 
2 Jo 

and the Wilson loop operator corresponds to the 
contour C enclosing the q, q paths from time 



to time T. It should include path-ordered par- 
allel transporters for the initial and final states, 
${x,x), $(y,y) with e.g. 

${x,x) = P ( dz^B^z). 

J X 

The calculation is simplified if choosing x = x, y = 
y which will be enough for our purposes here. To 
take into account the nonperturbative character of 
the interction it is convenient to work in the back- 
ground gauge formalism and write = b^ + 
where the a M represent the quantum fluctuations 
and b^ is a background field which is choosen such 
that the vacuum expectation value of the Wick 
ordered products of the vanish. Therefore, we 
may express the gluon correlator in terms of b^ 
only: 

(: G(x)G(y) :) - (: G b (x)G b (y) :), 
Gb,fi V = d^K - d v bn + g[bp, b v \. 

Expanding in powers of the background field b^ 
we may write the Wilson loop average as 

(W(C)}= J VaPe-fc dz » a » 

+ (l) j Va S c ^ J dz - P ^ (z ' ?J) ( 3 - 3 ) 
xb^{z)P^ a (z',z)b u (z') + ... 

= W + W 2 + ... 

and the transporter <P a is constructed with only 
the quantum field a. For the first term, Wq, the 
cluster expansion gives 

Wo = Z exp {<j>2 + higher orders) , 

= C F OL s r- x f dt + 0(v 2 ), 
Jo 

i.e., the Coulombic piece of the potential. (Z is 
a constant that, in particular, includes regulariza- 
tion). 

The evaluation of the first nontrivial piece, 
Wi is more complicated. It produces a correction 
to the Green's function, SG, which in the static 
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approximation is 

SG = -± J d 3 r J dV J nd0 J r\A0 

xG< S »(r(r),r)G< 8) (r,r')G< S »(r',r(0)). 

Here the G^ are the singlet, octet Coulombic 
Green's functions. We may then take matrix el- 
emets between Coulombic staes, \nl), and identify 
the ensuing energy shifts from the relation 

G = G { J } +5G ~ G { J\l-TSE nl ). 

T—*oo 

We then find the basic equation^ 18 ] , 

SE - = ™ / w / dm ' A{p) 

xY,(nl\r^ p( - - 1/2)r \k(8)) (3.4) 



E^-E n -p 



(k(8)\rle 



/„ip'(/3-l/2)r'| 



ll). 



The states |fc(8)) are eigenstates of the octet Hamil- 
tonian, with energy ; the E n are the Coulom- 
bic energies. Finally, A(p) is defined in terms of 
the correlators, being the Fourier transform of 

A(x) = D{x) + Dx{x) + x 2 d 2 D l {x)/dx 2 

and 

{g 2 :G ai {x)G Qj {Q) :) = £ 



5ijD(x) + XiX, 



d 2 D 1 



dx 2 



We may write, using Lorentz invariance, A{x) = 
f(x 2 /T 2 ), with T g the so-called correlation time. 
This will play an important role in what follows. 

We have now two regimes. If n. T = Tg 1 3> 
\E n \ the velocity tends to zero, and the nonlocal- 
ity also tends to zero as compared with the quark 
rotation period (which in the Coulombic approxi- 
mation would be 1/| E n \). We can now neglect, in 
Eq. (3.4), both E n , E\' as compared to po so we 
obtain 5E n i ~ (nl\U\nl) where 



/r roc -\ 

AdA / di/[-2D(A, i/)+£>i (A, !/)][• 

(3.5) 

At large r, and as this equation shows, we find 
U (r) ~ err. Here a can be related to T g and the 



gluon condensate. If e.g. we take an exponential 
ansatz for A(x), 

MT = -^^^0.32GeV. 

ov ^ (J 2 

For small r,! 16 ' 17 ! 

U(r) ~c + cir 2 . (3.6) 

This is different from the behaviour expected from 
the Leutwyler-Voloshin analysis which gives a be- 



haviour 



but one should understand that 



the present derivation holds for r — > but still 
T^ 1 ^> \E n \. It may be noted that the analy- 
sis based upon the potential U gives a very good 
description of heavy quarkonia states' 19 ' . 

We next get the matching between the 
two regimes! 18 ] p or this we now turn to the oppo- 
site situation, viz., T g ' 1 -C \E n \. Now we may ap- 
proximate A{x) ~ constant so that A(p) <~ &i(p) 
and Eq. (3.4) becomes 

7r(a s : G 2 :) 
SE ni = — — -<nZ|ri 



1 



if ( 8 ) - + MT 



|nZ), 
(3.7) 

which coincides exactly with the results of the 
Leutwyler-Voloshin analysis' 2,6 ! in the limit T g — > 
oo (fj,T — > 0). In fact, Eq. (3.7) allows us to esti- 
mate the finite size corrections to the NP effects, 
which improves still the agreement between theory 
and experiment^ 18 ! . 

4. RENORMALONS. SEMICLASSI- 
CAL UNDERSTANDING OF THE 
HEAVY QUARK POTENTIALS. 
SHORT DISTANCE LINEAR POTEN- 
TIAL AND SATURATION 
In the previous section we have shown how QCD 
can give a very satisfactory account of the heavy 
quarkonia spectra, particularly of the lowest ly- 
ing states; an understanding based on perturba- 
tive calculations supplemented by NP ones, in par- 
ticular those associated with the gluon conden- 
sate. Here we address two questions related to 
that. First, one may inquire about the connec- 
tion of renormalons with nonperturbative effects. 
Secondly, one can try to understand intuitively 
the potentials one finds. Finally, we will devote a 
few words to a speculation on a posible linear po- 
tential at short distances, and its connection with 
saturation. 
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FIGURE 4. One-gluon exchange, 
dressed with loops. 

Renormalons. Let us return to the one-gluon ex- 
change diagram, Fig. 1. If we dress the gluon prop- 
agator with loops (Fig. 4) then the corresponding 
potential, in momentum space, is 



V(k) = 



-AnC F 



Ait 



k 2 /3 log(fc 2 /yl 2 



(4.1) 



and we have substituted the one-loop expression 
for a s (k 2 ). The expression (4.1) is undefined for 
soft gluons, with k 2 ~ A 2 . As follows from the 
general theory of singular functions, the ambiguity 
is of the form cS(k 2 — A 2 ): upon Fourier transfor- 
mation this produces an ambiguity in the x-space 
potential of SV(r) — c[smAr]/r. At short dis- 
tances we may expand this in powers of r and find 



SV(r) ~ Co + C x r 2 + ... 



(4.2) 



The same result may be obtained with the more 
traditional method of Borel transforms' 20 ' 21 !. This 
coincides with the short distande behaviour of the 
nonperturbative potential U(r) as determined in 
refs. 13-17, and Eq. (3.6) here. [For applications to 
calculations of bound states, see ref. 22 and work 
quoted there]. 

The situation just described applies for 
states qq at short distances; but not so short that 
zero frequency gluons cannot separate the pair. If 
this last is the case, soft gluons do not resolve the 
qq pair and only see a dipole. The basic diagram 
is no more that of Fig. 1, but that of Fig. 5. The 
generated renormalon may then be seerJ 21 ' 23 ! to 
correspond to the contribution of the gluon con- 
densate in the Leutwyler-Voloshin mechanism. 




FIGURE 5. Emision and absorp- 
tion of a soft gluon collectively by 
a qq pair. 

Semiclassical picture We have seen that one can 
get a consistent QCD description of heavy quarko- 
nium ground states both for large and small T g . 
Here we will try and show how one can give an in- 
tuitive picture of what we have found I 21 ' . For this 
we consider a model for quarkonium to be that 
of a e + e~ pair inside a conducting cavity of ra- 
dius R ~ A~ x . The potential energy of the pair is 
given as an integral over space of the correspond- 
ing electric fields, 

V ( r ) = ^J d 3 r5!(r')£ 2 (r + r'). (4.3) 

In particular, the Coulomb potential is obtained 
when the £, correspond to point charges. If these 
fields are modified at long distances, this will give 
rise to a modification of this interaction also at 
small distances. In our case, the modification 
arises because, since the charges are confined, the 
integral in (4.3) should only be extended to r < R. 
Thus, 

which reproduces the quadratic term in (3.6). The 
constant term appears because now we cannot fix 
the Coulomb potential by requiring it to be zero 
at infinity. 

This calculation does not take into ac- 
count retardation effects. When these become im- 
portant, which is when the e + e~ pair is rotating 
very closely, the quadratic potential is wiped out 
and there remains a cubic one -again as in the 
QCD case. The situation is fully analogous to that 
of the ordinary Casimir effect! 24 ]. 
A linear potential at short distances? To finish 
this note we are going to speculate on the pos- 
sibility of a linear correction to the potential at 
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short distances. We have no proof of the exis- 
tence of such term, but we have three different 
indications for its existence. First of all we have 
the posibility that the QCD coupling saturates at 
long distances' 25 ] so that one has, 



a s (k 2 ) 



4-7T 



fc^o A)log[(fc 2 + M 2 )/A 2 



(4.4) 



with M ~ A (the possibility that M = A is sug- 
gested by the deep inelastic scattering evaluations 
of the second paper of ref. 25). This yields a linear 
potential correction when inserted in a Coulombic 
potential both at long and short distances. 

The second indication comes from lattice 
QCD calculations, where a linear correction to 
the short distance Coulombic potential is appar- 
ently seen' 26 ]. The third indication comes from 
the following intuitive argument' 21 '. Consider a 
simplified model according to which the chromo- 
electrostatic field of quarks is a correct zeroth- 
order aproximation only so far as it exceeds some 
critical value of order A 2 : £ 2 > A 4 , while weaker 
fields do not penetrate the vacuum because of its 
specific, confining properties. From this condi- 
tion we get an estimate of distances R CY where 
the chromo-electrostatic field of quarks is strongly 
modified: 



<A 4 



where for simplicity we have neglected the effect 
of the running of a s (r~ 2 ). 

The corresponding change in the potential 
is then of order 



SV 



t\l 2 rA 2 



(4.5) 



i.e., we get a leading correction linear in r to the 
potential at short distances. 

It is not easy to see how one could get 
a handle on this linear potential. The agreement 
between the orthodox QCD calculations and ex- 
periment is so good (see above and e.g. refs. 6, 
8, 10) that there seems to be little room for (4.5). 
The saturation modification of a s would also be 
masked by the errors in A. Perhaps lattice calcu- 
lations may give a hint, as they seem to be doing 
already' 26 ! . 
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